Intermittent dislocation flow in viscoplastic deformation 
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The viscoplastic deformation (creep) of crystalline materials under constant stress involves 
the motion of a large number of interacting dislocations |jl|]. Analytical methods and sophisti- 
cated 'dislocation-dynamics' simulations have proved very effective in the study of dislocation 
patterning, and have led to macroscopic constitutive laws of plastic deformation [^]^. Yet, 
a statistical analysis of the dynamics of an assembly of interacting dislocations has not hith- 
erto been performed. Here we report acoustic emission measurements on stressed ice single 
crystals, the results of which indicate that dislocations move in a scale-free intermittent fash- 
ion. This result is confirmed by numerical simulations of a model of interacting dislocations 
that successfully reproduces the main features of the experiment. We find that dislocations 
generate a slowly evolving configuration landscape which coexists with rapid collective rear- 
rangements. These rearrangements involve a comparatively small fraction of the dislocations 
and lead to an intermittent behavior of the net plastic response. This basic dynamical pic- 
ture appears to be a generic feature in the deformation of many other materials p!o|-p!^. 
Moreover, it should provide a framework for discussing fundamental aspects of plasticity, 
that goes beyond standard mean-field approaches that see plastic deformation as a smooth 
laminar flow. 



Whenever dislocation glide is the dominant plas- 
tic deformation mechanism in a crystalline material, 
we observe a constant strain-rate regime usually de- 
scribed by Orowan's relation 7 ^ Pmbv. Here, the 
plastic strain-rate of the material 7 is simply re- 
lated to average quantities such as pm, the density 
of mobile dislocations, and w, their average veloc- 
ity along the slip direction (parallel to the Burgers 
vector b) Q. Transmission electron micrographs 
of plastically deformed materials display, on the 
other hand, complex features such as cellular struc- 
tures and fractal patterns H,^, which are the fin- 
gerprint of a complex multiscale dynamics not ap- 
propriately accounted for by the mean-field char- 
acter of Orowan's relation. In addition, rapid slip 
events ||l^ have been observed in the plastic de- 
formation of various metals and alloys ||Tl|,|l^, and 
in the Portevin-LeChatelier effect |l4|. We believe 
that formulating plastic deformation as a nonequi- 
librium statistical mechanics problem |l4[ requires 
a substantial understanding of basic collective dis- 
location dynamics. 

Experimentally, the complex character of collec- 
tive dislocation dynamics can be revealed by acous- 
tic emission measurements. The acoustic waves 
recorded in a piezoelectric transducer disclose the 
pulse-like changes of the local displacements in 
the material during plastic deformation, whereas 
a smooth plastic flow would not be detected [|15|. 



Thus, this method is particularly useful for inspect- 
ing possible fluctuations in the dislocation velocities 
and densities. 

Ice single crystals can be used as a model mate- 
rial to study glide dislocation dynamics by acoustic 
emission |l^,0 due to the following reasons: (i) 
Transparency allows direct verification that acous- 
tic emission activity is not related to microcrack- 
ing. (ii) Within the range of temperature and stress 
corresponding to our experimental conditions, dif- 
fusional creep is not a significant mechanism of in- 
elastic deformation [p^ which, in hexagonal ice sin- 
gle crystals, occurs essentially by dislocation ghde 
on the basal planes along a preferred slip direc- 
tion, (iii) An excellent coupling between sample 
and transducer can be obtained by fusion/freezing. 

Uniaxial compression creep experiments are per- 
formed on artificial ice single crystals, employing 
several steps of constant applied stress. We observe 
an intense acoustic activity, exhibiting a strong in- 
termittent character (see the inset of Figure |l|). 
We measure the energy associated to each acous- 
tic burst and analyze its statistical properties. In 
Fig. we show that the probability distribution of 
energy burst intensities exhibits a power law be- 
havior spanning several decades. This fact is an 
indication that a large number of dislocations move 
cooperatively in an intermittent fashion. A sim- 
ilar behavior has been observed in the Portevin- 



LeChatelier effect ||l3|, a plastic instability where 
the intermittent flow is ruled by the interaction of 
dislocations and diffusing solute atoms. The inter- 
mittency observed in our case is different as we only 
have interacting dislocations subject to an external 
stress, without any other alien element interfering 
with their dynamics. 
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FIG. 1. Statistical properties of the acoustic energy 
bursts recorded in ice single crystals under constant 
stress. The main figure shows tlie distribution of en- 
ergy bursts for tlie different loading steps. Tfie cutoffs 
depend slightly on the applied stress, but the power law 
exponent remains the same. The fit yields an exponent 
te = 1.60 ± 0.05. Inset, a typical recorded acoustic sig- 
nal. The creep experiment is performed at T = — 10°C 
under different constant compression stresses ranging 
from a = 0.58MPa to a = 1.64MPa. The stress is 
apphed at an angle with respect to the c-axis, giving 
rise to a small resolved shear stress acting across the 
ghde plane (ae = 0.030MPa-0.086MPa). The frequency 
bandwidth of the transducer was 10 to lOOkHz. The 
amplitude range between the minimum and the maxi- 
mum recordable thresholds is 70dB, that is, a range of 
seven orders of magnitude in energy. 

Multiscale properties and pattern formation are 
ubiquitous in plastic materials and we expect that 
the large dynamical fluctuations observed in ice 
compression experiments are also a prevalent fea- 
ture of plastic deformation micromechanics. This 
general picture implies that the experimental phe- 
nomenology should be reproducible in simple nu- 
merical simulations of discrete dislocation dynam- 
ics that preserve the relevant characteristics of the 
system under study. As in previous dislocation dy- 
namics models [H-M , we consider a two-dimensional 
cross-section of the crystal (that is, the xy plane), 
and randomly place N straight-edge dislocations 
gliding along a single slip direction parallel to their 
respective Burgers' vectors b (that is, the x direc- 



tion). This implies that we have point-like dislo- 
cations moving along fixed lines parallel to the x 
axis. This simplification effectively describes ma- 
terials like ice crystals that, owing to their strong 
plastic anisotropy, deform by glide on a single slip 
system. An edge dislocation with Burgers vector b 
located at the origin gives rise to a shear stress Us 
at a point r = (x, y) of the form 



27r(l - z/)(a;2 -|-y2) 
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where /i is the shear modulus and u is the Poisson 
ratio jl]. This long-range stress is responsible for 
mutual interactions among dislocations which are 
important in all dislocation dynamics models . 
Under a constant external stress cTg, dislocation i 
performs an overdamped motion along the x direc- 
tion described by the following equation 



(2) 



where "q is the effective friction and hi is the Burg- 
ers vector. Throughout the simulations, we will 
be dealing with dimensionless quantities after set- 
ting the distance scale 6=1, and the time scale 
to = ilb/{ti/2T:{l-v)) = l. 

Other essential ingredients commonly introduced 
in most dislocation dynamics models iQ-Q are the 
mechanisms for dislocation (i) annihilation and (ii) 
multiplication, (i) When the distance between two 
dislocations is of the order of a few Burgers vec- 
tors, the high stress and strain conditions close to 
the dislocation core invalidate the results obtained 
from a linear elasticity theory (equation (0)). In 
these instances, phenomenological nonlinear reac- 
tions describe more accurately the real behavior of 
dislocations in a crystal. In our model, we annihi- 
late two dislocations with opposite Burgers vectors 
when the distance between them is shorter than 25. 
(ii) The activation of Frank-Read sources |^ is ac- 
cepted as the most relevant dislocation multiplica- 
tion mechanism under creep deformation. The ac- 
tivation of Frank-Read sources have been observed 
in ice crystals, along with more specific multipli- 
cation processes @. Because an accurate multi- 
plication mechanism cannot be simulated in a two- 
dimensional point-dislocation model, we have im- 
plemented various phenomenological procedures. A 
simple way of taking into account the new disloca- 
tion loops generated at different sources within the 
crystal is the random introduction of opposite sign 
dislocation pairs in the cross-section under consid- 
eration. The rate and frequency of this creation 
process depend solely on the external stress Ce. 
We have checked that other multiplication rules in 
which the creation rate depends on the local stress 
yield similar results. 
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FIG. 2. Snapshot of the total stress field and dislocations arrangement m a numerical simulation with Va — 0.025. 
We observe metastable structure formation (dipoles and walls) and the associated stress field which goes from light 
blue (low values) to dark blue (high values). The complex low-stress pathways joining dislocations are the result of 
the anisotropic elastic interactions. Dislocations moving at low velocities (lower than Vcr) are depicted in magenta, 
while those moving at higher velocities are depicted in black. Most dislocations in walls are moving slowly, whereas 
isolated ones tend to move at higher speed. In the simulations, an initial number of A^o = 1500 dislocations are 
randomly distributed on a square cell of size L = 300. Their Burgers' vectors are randomly chosen to be +b or —b 
with equal probability. In the absence of external stress, we first let the system relax until it reaches a metastable 
arrangement. The number of remaining dislocations is then N ~ 700. At this point, we apply a constant shear 
stress and study their evolution. To avoid the discontinuities arising from truncating long-range elastic interactions 
in Eq. (^), we resort to the Ewald summation method. We have thus exactly taken into account the interaction of 
a dislocation with all the infinite periodic images of another dislocation in the same cell. 
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FIG. 3. Statistical properties of dislocation veloci- 
ties and density obtained in numerical simulations, a) 
Probability density of the individual dislocation veloc- 
ities \vi\ in the stationary state. The black line is a 
least-squares fit of the intermediate data points. We 
obtain an exponent u = 2.5 for the scaling of the in- 
termediate velocities, b) Time evolution of the total 
number of dislocations in the cell, c) Fraction of fast 
moving dislocations for a given run of the simulations. 

In the initial stage of the simulations, the sys- 
tem relaxes slowly and the average strain-rate de- 
creases until it reaches a constant value which de- 
pends on the applied stress, that is the stationary 
creep regime. By monitoring the activity in a single 
run, we observe that most dislocations are arranged 
into metastable structures (such as walls and cells) 
moving at a very slow rate (see Figure ^ and the 
Supplementary Information). A smaller fraction of 
dislocations, however, move intermittently at much 
higher velocities, giving rise to sudden increases of 
the plastic strain. In Figure^), we plot the prob- 
ability density of finding an individual dislocation 
moving with a velocity \vi\. Already at the level 
of individual dislocations, the velocity distribution 
P(|wi|) is quite broad and exhibits power law be- 
havior for velocity values larger than the external 
stress induced velocity = bae/rj (we have consid- 
ered two values of Va ■ 0.0125 and 0.025). 

In the presence of a large number of dislocations, 
the acoustic signal that is detected experimentally 
is due to the superposition of the waves emitted 
by each moving dislocation In particular, it 



has been shown that a single dislocation perform- 
ing a sudden movement with velocity vq for a short 
time r gives rise to an acoustic wave whose am- 
plitude is proportional to vq 15 |. The high am- 



plitude pulses detected in the experiment cannot 
be ascribed to uncorrelated emissions from each in- 
dividual dislocation, but rather to the cooperative 
motion of several dislocations occurring, for exam- 
ple, after the activation of a Frank-Read source, or 
if a dislocation cluster breaks apart. To gain fur- 
ther insight into this behavior, in Figures yb) and 
^), we show the evolution of the total number of 
dislocations as well as that of the fast moving dis- 
locations, that is, dislocations moving faster than if 
they were moving under the only action of the ex- 
ternal stress, Vi > (other threshold values yield 
equivalent results). After the injection of a few new 
dislocations or the annihilation of a dislocation pair, 
several other dislocations start to move and rear- 
range, not necessarily in the close vicinity of the 
triggering event. To quantify this effect, we mea- 
sure the collective velocity V = J2 oi the fast 
moving dislocations and define the acoustic energy 
as £' = V^^ [0. In the inset of Figure^, we see that 
the signal E{t) consists of a succession of intermit- 
tent and pronounced bursts, each one signalling the 
onset of collective dislocation rearrangements. The 
slow and continuum motion of dislocation struc- 
tures (vi < Va) is not considered as it only sets up a 
background noise signal which cannot be detected 
experimentally. In Fig. U) we show that the distri- 
bution of energy bursts, obtained by sampling the 
signal over different times and realizations, decays 
very slowly, spanning various decades. For interme- 
diate values, the distribution shows an algebraic de- 
cay with an exponent te = 1.8 ± 0.2, in reasonable 
agreement with experiments (see Fig. |l|). The max- 
imum number of dislocations we can handle in our 
simulations severely restricts the maximum value of 
the signal in the system. Consequently, the exten- 
sion of the power law regime grows with the num- 
ber of dislocations, but is eventually bounded by 
the different nature of the extremes statistics. 

The broadly distributed intermittency resulting 
from the statistical analysis can be interpreted as a 
nonequilibrium transport phenomenon with scaling 
properties. Scaling behavior is usually associated 
with the proximity of a critical point and is charac- 
terized by a high degree of universality. Critical ex- 
ponents only depend on a few fundamental parame- 
ters such as the effective dimensionality, and the ba- 
sic symmetries of the system. Specifically, the rea- 
sonable quantitative agreement between our model 
and the experimental data is due to the strong plas- 
tic anisotropy present in ice single crystals, that 
can thus be well described by a two dimensional 
model. Although the quantitative results we ob- 
tain should be restricted to the case of single slip 
systems, the general features of the observed inter- 
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mittent flow regime appear to be generic in plas- 
tic deformation [ pl]JT2| . The particular miiversal- 
ity class will depend on the effective dimensionality 
of each specific material as well as the particular 
deformation mechanism, and one can still expect 
to obtain the right critical exponents from simpli- 
fied models if the relevant symmetries are correctly 
taken into account. 




FIG. 4. Statistical properties of the energy bursts ob- 
tained in numerical simulations. In the model, we de- 
fine the acoustic energy E = from the collective 

velocity V = X/i I""*! '^^ f"'^^ moving dislocations. 
The time evolution of E for a given simulation run is 
plotted in the inset. The energy distribution averaged 
over time and over a hundred realizations with differ- 
ent initial conditions is plotted in the main graph. The 
black line represents the best fit of the intermediate data 
points. We obtain an exponent te = f .8 ± 0.2 for both 
values of the external stress. 

The close-to-criticality nonequilibrium frame- 
work motivates possible analogies with elastic man- 
ifolds driven in random media such as fluid flow 
in porous media, vortices in superconductors, and 
charge density waves |^^. In all of these systems, 
a critical value of the driving force separates a 
static regime from a moving one, and scaling is 
observed only close to this point. The intermit- 
tent behavior close to criticality is in these cases 
associated to static random heterogeneities which 
exert space dependent pinning forces on the mov- 
ing objects. On the contrary, the present disloca- 
tion dynamics model, as well as the ice experiment, 
does not contain any quenched source of pinning 
forces. Dislocation themselves, through the vari- 
ous structures such as dipoles and walls, are self- 
generating a pinning force landscape in which the 
dynamics is virtually frozen; i.e. a slow dynamics 
state. Creation and annihilation mechanisms, often 
triggered by the presence of unsettled dislocations. 



allow the system to jump between slow dynamics 
states through bursts of activity. This behavior is 
reminiscent of driven-dissipative self-organizing sys- 
tems that achieve criticality in the limit of a very 
slow driving [pT| . 

The emerging scenario poses many basic theoret- 
ical questions. Is there a critical stress value be- 
low which the system decays into a slow dynam- 
ics state? Can the large dynamical fluctuations be 
associated with a diverging response function and 
thus, be related to a critical point? Can we derive 
scaling laws relating the various observed exponents 
as suggested by the theory of critical phenomena? 
The answers to all these questions pave the way to 
the nonequilibrium statistical theory of dislocation 
motion that is needed for a deeper understanding 
of the micromechanics of plastic deformation. 
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